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Abstract 

Recently, it has been suggested that large- N c orbifold equivalences may be applicable to cer- 
tain theories with chemical potentials, including QCD, in certain portions of their phase diagram. 
When valid, such an equivalence offers the possibility of relating large- N c QCD at non-zero baryon 
chemical potential, a theory with a complex fermion determinant, to a related theory whose 
fermion determinant is real and positive. In this paper, we provide a test of this large N c equiva- 
lence using a holographic realization of a supersymmetric theory with baryon chemical potential 
and a related theory with isospin chemical potential. We show that the two strongly-coupled, 
large- N c theories are equivalent in a large region of the phase diagram. 
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1 Introduction 

Understanding QCD at non-zero baryon density is an important goal, both for its intrinsic 
interest and for applications such as the structure of neutron stars and the mechanism of core- 
collapse supernova. Due to the notorious sign problem^ we lack generally effective methods for 
performing numerical simulations of gauge theories with a baryon chemical potential. When a 
non-zero baryon number chemical potential is present, the determinant of the Euclidean Dirac 
operator is no longer positive and standard Markov-chain Monte-Carlo methods are not applicable. 
Although many schemes have been proposed to address the sign problem [U GJ EJ [H El [7] , it is 
fair to say that no fully satisfactory solution has been found. At the same time, condensed matter 
phases of several QCD-like theories which do not suffer from the sign problem have been studied 
numerically, in the hope that one may extract lessons about strongly interacting finite density 
systems which will also apply to QCD at finite baryon density. Examples include SU(2) Yang- 
Mills (YM) with even numbers of fundamental flavors SU(N C ) YM with adjoint fermions 
[9], and QCD with an isospin chemical potential [10|. Illj . However, there is no solid argument 
delineating the extent to which these theories can reproduce properties of QCD with a baryon 
chemical potential. 

In recent years, it has been understood that a network of large- N c equivalences relate various 
non- Abelian gauge theories with differing gauge groups and matter content [121 El O EG3 IS] • 
These equivalences, which are generated by appropriate orbifold projections, relate the leading 
large N c behavior of connected correlators of specific classes of observables. The large- N c equiv- 
alences are valid provided certain symmetry realizations are satisfied |17j . For example, SU(N C ) 

5 More properly, this should be called a phase problem. See Appendix \A\ for a brief summary. 
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0{2N C ) 
Nf fundamentals 



chemical potential \xp 



U(N C ) 
Nf fundamentals 
baryon chemical potential [Ib 



U(N C ) 
Nf fundamentals 
isospin chemical potential /i/ 



Figure 1: An orbifold projection acting on a parent 0(2N C ) Yang-Mills theory with JVf fundamental 
representation Dirac fermions and a flavor-singlet chemical potential fj,p may generate a U(N C ) 
daughter theory with Nf fermions and a baryon chemical potential fiB (left) or, provided Nf is 
even, the same U(N C ) theory with an isospin chemical potential [ii (right) [21J. In the parent 
0(2N C ) theory with even Nf, there is no distinction between a baryon or isospin chemical potential. 



and SO(2N c ) Yang-Mills theories have coinciding large N c limits of all Wilson loop expectation 
values (as well as connected correlators of real parts of Wilson loops), provided charge conjugation 
symmetry is not spontaneously broken in the SU(N C ) theory |18) . 

We will be concerned with QCD-like theories containing fundamental representation fermions 
and non-zero chemical potentials. Specifically, we will discuss: 

1. SO(2N c ) Yang-Mills with Nf Dirac fundamental representation fermions and a non-zero 
fermion chemical potential /ip, under which all Nf flavors have charge +1. For brevity, we 
will denote this theory as SO(2N c )p. 

2. SU(N C ) Yang-Mills with Nf fundamental representation fermions and a non-zero baryon 
chemical potential ^b, under which all Nf fermion flavors have charge +1. For brevity, we 
will denote this theory as QCD#. 

3. SU(N C ) Yang-Mills with Nf fundamental representation fermions, with Nf even, and a non- 
zero isospin chemical potential \xi, under which half the fermion flavors have charge +1 and 
half have charge —1. For brevity, we will denote this theory as QCD/. 

Although QCD# suffers from a sign problem, this is not the case for either QCD/ or SO(2N c )p, 
as both of these theories have a real and positive fermion determinant ¥U\ . 

As figure [T] schematically depicts, starting from the SO(2N c )p theory one choice of orbifold 
projection yields QCD#, while a different choice yields QCD/j^] Based on this observation, it has 
recently been suggested that large N c equivalences may relate suitable observables in the parent 

6 Strictly speaking, the orbifold projection maps a parent theory with SO(2N c ) gauge group to a daughter U(N C ) 
gauge theory. But the difference between U(N C ) and SU(N C ) theories is sub-dominant in the large N c limit. Note 
that in a U(N C ) theory with chemical potential, the U(l) part of the gauge field is to be fixed at infinity. In finite 
volume, the theory should be defined with Dirichlet boundary conditions on the U(l) gauge field, not periodic. (In 
practice, for lattice simulations, it is more convenient to simply use the SU(N C ) theory.) 
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S0(2N c )p theory to corresponding observables in either QCD# or QCD/ [2T1 [22J, [23] . In portions 
of the phase diagram where both equivalences are valid (if such regions exist), this implies that 
one may obtain quantitative information about large- N c QCD with a baryon chemical potential 
from studies of the same theory with an isospin chemical potential, thereby circumventing the 
sign problemj^] 

When N c — > oo with N{ fixed, a comparison of planar Feynman diagrams in the parent 
SO(2N c ) and daughter SU(N C ) theories shows that they coincide [H 13 |2H [22] . (For other ap- 
proaches see [E] •) This analysis is unaffected by the presence of a non-zero chemical potential. Co- 
inciding perturbative expansions does not, however, necessarily imply a genuine non-perturbative 
equivalence. Necessary conditions for a valid equivalence include a requirement that the symme- 
tries used to define the orbifold projection not be spontaneously broken in the parent theory [15j . 
Since the projection leading to QCD b is generated by a combination of a gauge transformation and 
a U(1)f phase rotation, this projection can only lead to a valid large- N c equivalence in portions 
of the phase diagram where the U{1)f global symmetry associated with net fermion number is 
unbroken. In other words, a large N c equivalence relating QCD with baryon and isospin chemical 
potentials can only apply to portions of the phase diagram in which fermions do not condense to 
form a superfluid. In simpler examples, analogous conditions on symmetry realizations are both 
necessary and sufficient conditions for the validity of large iV c equivalences [15]; whether this is 
the case in the present example is not yet clear. 

In this paper, we use gauge/gravity duality to test the validity of analogous possible large iV c 
equivalences relating supersymmetric generalizations of the above theories. Although the trio of 
theories we consider will not include QCD itself, the arguments of refs. [2TJ [221 E3] are equally 
applicable to the supersymmetric theories we consider. By considering supersymmetric theories, 
and using holographic methods, it will be possible to examine relations between theories with 
different types of chemical potential directly in the limit of strong coupling (and large N c ) using 
simple analytic methods. We will find that, in a large region of the phase diagram with no 
spontaneous breaking of flavor symmetries, large N c equivalences between our theories are valid. 

Our holographic construction involves an orbifold and an orientifold projection of the D3/D7- 
system, with N c D3 branes, iVf D7 branes, and always Nf <C N c . At large- iV c and strong 't Hooft 
coupling, the low energy theory on the D3 branes is described by classical type IIB supergravity 
in AdSz, x S* 5 , with probe D7 branes wrapping an S 3 in the S 5 [24J. The projections act on the 
geometry, changing it to AdS§ x RP 5 . An isospin chemical potential in the original theory is 
described by a particular configuration of the gauge field on the D7 brane. We show that after the 
projection the isospin chemical potential becomes a baryon chemical potential. Figure [2] illustrates 
the connections between the corresponding field theories. We prove that, provided the projection 
symmetries are not broken, the equations of motion of the D7 branes coincide in both theories. 
As we discuss, this implies that the conjectured large N c equivalences are valid in these theories 
in those regions of the phase diagram where flavor symmetries are unbroken and no additional 
fields become active. 

The paper is organized as follows: in §[2] we review in more detail the proposed equivalences 
involving large- N c QCD. In §|3]we provide a holographic realization for our supersymmetric gener- 
alization. Finally, in §[4] we discuss the regime of validity of the equivalence, and the consequences 
for the phase diagram. 

7 Readers should refer to section 2, and refs. |21l 1221 [23] . for more details and more nuanced discussion. 
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SU(2N C ), 2N{ fundamentals, 3 adjoints 
isospin chemical potential [ij 



SO(2N c ), N{ fundamentals, 1 adjoint, 2 symmetries 
chemical potential \ip 

SU(N C ), Nf fundamentals, 1 adjoint, 2 symmetries 
baryon chemical potential \xb 

Figure 2: M = 1 supersymmetric theories related by orbifold projections. In the SO(2N c ) theory, 
there is no distinction between a baryon or isospin chemical potential. 

2 Orbifold projections for large-A^ c QCD 

Consider an SO(2N c ) Yang-Mills theory coupled to Nf fundamental representation Dirac 
fermions. The Lagrange density is 

1 Nf 

£so = -r-r- tr F lu +^2*Pa (7 M ^ + m g + fi F 7°) i> a , (1) 
^ 9 so ' tt=1 

where is the field strength of the SO(2N c ) gauge field A^, = d^ + A^, ip a is a Dirac fermion 
in the vector representation of SO(2N c ), and m q and fip are the quark mass and fermion chemical 
potential, respectively. Because the gauge field is real, the Dirac operator V = {^^D^+mq+^p 7 ) 
satisfies (C^)D (Oys) -1 = D*, where C is the charge conjugation matrix defined by Cj^C^ 1 = 
— = — 7^0 If v is an eigenvector of the Dirac operator T> with an eigenvalue A, T>v = Xv, then 
(C75) _1 f* is another eigenvector of T> with eigenvalue A*, and is linearly independent of v even 
when A is real. (See § 2.3 of ref. [20].) Therefore the determinant of D is always real and positive, 
implying that standard Markov chain Monte-Carlo simulation techniques may be used [21] . 

When m q = fip = 0, the Lagrangian ([!]) has a manifest SU(Nf)L x SU(N{)r x U(l)p x 
U(1)a flavor symmetry, just like SU(N C ) QCD. However, the flavor symmetry of the theory 
is actually larger than this due to the fact that SO(2N c ) is a real gauge group; classically it 
extends to U(2Nf) [25, 26J. The axial U(1)a C U(2Nf) is anomalous, and at the quantum 
level the (continuous part of the) flavor symmetry is SU(2Nf), which spontaneously breaks to 
SO(2Nf) ^ SU(Nf)v due to the formation of a chiral condensate ('tpip). The resulting massless 
Nambu-Goldstone bosons span the SU(2Nf) / SO(2Nf) coset space. In contrast to QCD, some 
of these Nambu-Goldstone bosons, which we will refer to as baryonic pions, are charged under 
U(l)p. Ordinary pions are created by operators that look like i^a'Jbipb, while baryonic pions are 
created by color-singlet operators of the form t/j^C^tptj and 'ipaCj^f . 

8 We use (++++) metric signature and Hermitian gamma matrices. 
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2.1 From SO(2N c ) F to QCD with a baryon chemical potential 

To perform an orbifold projection, one identifies a discrete subgroup of the symmetry group 
of the parent theory, which for us is the SO(2N c )p theory, and then removes all of the degrees of 
freedom in the parent theory which are not invariant under the chosen discrete symmetry. This 
yields a daughter theory, which will turn out to be large- N c QCD. 

The required orbifold projection is a Z2 subgroup of the SO(2N c ) gauge x U{1)f flavor sym- 
metry of the SO(2N c ) p theory. To define the orbifold projection, which we will denote as Vb, take 
<^27V C 6 SO(2N c ) to be given by J2N C = i&2 <8> 1jv c , where lj\r denotes an N x N identity matrix. 
The group element J2N C generates a Z4 subgroup of SO(2N c ). Next, let u = e i7T ^ 2 G U(1) F denote 
the phase which generates a Z4 subgroup of U{\)f- The discrete symmetry which will define the 
orbifold projection acts on on the fields A^, ip a as 

A-p, -> J2N c \J2N c i V>a ->■ UJ-ZN^O, ■ (2) 

Since J$ N = —l2N c , and (luJ2N c ) 2 = +l27V c! this symmetry transformation generates a Z2 sub- 
group of SO{2N c ) x U(1) F . 

The action of the orbifold projection on the basic fields is 

VbA^ = \ (A„ + J 2 N C A^2N C ) . Vb^ = \U + J2N c ^K n ]) (3) 



where, for later convenience, we have defined a matrix K^ 1 = ilN { acting on flavor indices. To 
display the action of the projection more explicitly, it is convenient to block-decompose the gauge 
and fermion fields. The gauge field A^ may be written in terms of four N c x N c blocks as 



where fields marked with an 'A' or i S > superscript are anti-symmetric or symmetric matrices, 

1^, and are even while B^, and C£ 



respectively. Under the Z2 symmetry transformation (|3|), A^ , and are even while B^, and 



are odd, so the orbifold projection sets B^ = = 0. Hence 



A A -n s 



If one defines the unitary matrix 



/> 1 ( l v ' ''' v ' s (0) 



then 



A u 



PVbA.P- 1 = ( "J* ) , (7) 



where = A A + iD^ is a U(N C ) gauge field. At large N c , we can neglect the difference between 
U(N C ) and SU(N C ) up to l/iV c 2 corrections. 
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We can split the 2A r c -component fundamental fermions of the SO(2N c ) theory into two N c - 
component fields, 

*-($)• (8> 

and then we use the matrix ([6]) to change basis. This yields 

(9) 

where ip± = (ipi ± 1^)2) I \f2. From eq. ([7]), one sees that V+ an d ip- transform as fundamental 
and antifundamental representations under SU(N C ), respectively. After the projection, only ip + 
survives. 

If we take the Lagrangian of the parent theory and apply the orbifold projection, it becomes 

, iVf 

£ = 7-2- Tr J% + V A a + m, + ^ B7 4 ) A a , (10) 

where J 7 ^ is the field strength of the SU(N C ) gauge field A^ = A^ + iD^, = + A^, 
X a = \/2i{)+, and the gauge coupling is given by g 2 su = g 2 so . 

(p) 

In the large- A^ limit for fixed iVf, connected correlation functions of operators 0\ in the 
parent SO theory which are invariant under the projection symmetry, and their counterparts 
in the daughter SU theory which are formed from the projected fields, coincide to all orders in 
perturbation theory [12], 

(o?o?...) p = (o?o? ..->* (11) 

The baryonic pion fields do not survive the projection, so there is no equivalent to them in the 
daughter theory. 

2.2 From SO(2N c ) F to QCD with an isospin chemical potential 

When the number of flavors in the parent SO(2N c ) theory is even, Nf = 2k, it is also possible to 
define a projection which yields large- N c QCD with an isospin chemical potential. The projection 
for the gauge field is the same as in eq. Q, but we now choose a different orbifold action on the 
flavor indices of the fermions. Let us write the fermions using N c x A/f-component fields as 

*-($$)■ 

In this basis, the orbifold action is 

f/>-> Jzn tI>J£- (13) 
This transformation also generates a 7L<i group. The action of the orbifold projection Vi is 

PjA^ = \ (Ap + J 2Nc A^ J~l) , Vii> = § + J2 N JJ 2 - k 1 ) ■ (14) 
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Defining tp± = (V>± =F iij)±)/\/2 and £± = (i/j±^ ± iij)±/)/y/2, one sees that <p± survive while £± 



is eliminated by the projection (14). Since tp + and cp- couple to Aa and A?, respectively, the 



fermionic part of the action of the daughter theory can be written as 

k 

£ £ A? ( 7 ^ M + m ± „ 7 4 ) A?, (15) 
/=i ± 

where A:jf = \/2tp]l , A_ = \/2 (ip^) c , and we have now written the flavor index (/) = 1, • • • , fc 
explicitly. This theory has an isospin chemical potential [ii = 2fi. 

2.3 Validity of large- iV c equivalences and their application to the sign problem 

The perturbative proof of the parent-daughter equivalence with isospin chemical potential 
is valid also when quark loops are included in planar diagrams, so it is possible to extend the 
analysis to include Nf/N c corrections. However, this is not possible for the projection to a theory 
with baryon chemical potential. The difference stems from the properties of the projection in 
the flavor sector, while the projection to isospin chemical potential is performed using a regular 
representation [121 [13] 

trJ 2fc = 0, Jf fc = ±l 2fe , (16) 

these conditions are not satisfied for the representation used to do the projection to the theory with 
a baryon chemical potential, where we have used -f^fc instead of J2k- Only diagrams containing a 
single quark loop produce the same result in parent and daughter theories. 

To go beyond the perturbative proof of the equivalence one needs to do a careful analysis of the 
necessary and sufficient conditions that must be obeyed for the equivalence to hold. A necessary 
condition is that the projection symmetry not be spontaneously broken in the parent |15j . The 
U(l) b symmetry, which is used for the projection from the SO(2N c ) theory to QCD with a baryon 
chemical potential, breaks to Z 2 when the baryonic pion condenses (e.g., when \x > m 7r /2 at zero 
temperature). Therefore, the parent-daughter equivalence can hold only at smaller values of the 
chemical potential^] On the other hand, the projection symmetry to obtain QCD with isospin 
chemical potential should not be spontaneously broken for any fj,; in this case condensation of 
baryonic pions in the parent theory is mapped to pion condensation in the daughter theory. 

Clearly, if it were possible to show that these equivalences hold nonperturbatively, they would 
be very useful because one would be able to derive properties of a large-iV c QCD theory with 
baryonic chemical potential from a SO(2N c ) theory or from large- N c QCD with isospin chemical 
potential, both of which are free of the sign problem. This could also explain why the phase 
quenching approximation in QCD is quite good — for a certain class of operators (e.g., the chiral 
condensate), the phase quenching approximation becomes exact in the large- iV c limit The phase 
quenching approximation for the chiral condensate is exact in the chiral random matrix model 
[39, 22]. The orbifold equivalence, if true, would ensure that the phase quenching approximation 
in QCD is exact for a large class of observables in the large- N c limit, even beyond the parameter 
region where the chiral random matrix model is valid (the "e-regime"). 



9 Note that the chemical potential at which baryonic pions condense is temperature dependent, and should 
increase with increasing temperature. 

10 Note that (for Nf even) dropping the phase of the fermion determinant turns the functional integral for QCDs 
into that for QCD/. 
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To provide a nonperturbative proof of the orbifold equivalence in QCD with chemical potentials 
is beyond the scope of this paper. However, in the following section we will show that analogous 
equivalences hold in a class of supersymmetric cousins of QCD which have gravity duals. 

3 A holographic realization 

It is possible to build a simple supersymmetric model where an isospin chemical potential 
is projected into a baryon chemical potential. The model is one of the examples mentioned in 
ref. |27| . based on the description of N = 2 theories from D4 branes suspended between NS5 
branes [28]. Flavor can be introduced by adding D6 branes. We will start with a configuration 
whose low energy limit on the T-dual D3 branes is N = 4 U(2N C ) super Yang-Mills plus 2iVf 
hypermultiplets in the fundamental representation, so the flavor group is U(2N{). In the T-dual 
configuration the flavor branes are D7's, and we will work in the 't Hooft limit of Nf/N c <C 1 so 
we can neglect their backreaction just as in the D3/D7 system of ref. [23]. We then introduce an 
orientifold plane to produce an SO(2N c ) theory with USp(2N{) flavor group and then finally do 
a Z2 orbifold projection that reduces it to U(N C ) with U(N{) flavor group. We will show that 
an isospin chemical potential in the original U{2N C ) theory is projected to a baryon chemical 
potential in the U(N C ) theory and discuss when the two theories are equivalent. 

3.1 Orientifold and orbifold projections 

The construction in type IIA theory consists on a set of 2N C D4 branes wrapping a circle in 
the x 6 direction and intersecting two 06 + planes at opposite sides of the circle. In addition, there 
is a NS5 brane at each orientifold point and 2N{ D6 branes parallel to the 06 planes: 
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Since the 06 planes are positively charged, Ramond-Ramond (RR) tadpoles do not cancel and 

the f3 function for the 't Hooft coupling is positive. However, in the 't Hooft limit Nf <C iV c , the 

/3 function is suppressed by Nf/N c at large N c . So to leading order in Nf/N c we can neglect the 

tadpoles and consider the _D6's and 06's as probes. 

This brane setup has as a T-dual a configuration involving D3 and D7 branes. The two 06 

planes map to a single 07 plane and the NS5 brane to a Z2 singularity localized at x 6 = x 7 = 
x 8 _ x 9 _ o: 
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The geometric effect of the Z2 action is a reflection in the transverse directions. The orientifold 
projection 0' = QR^(— 1) Fl involves worldsheet parity reversal f2, a reflection R45 in the x 4 and 
x 5 coordinates, and (— 1) Fl acts as —1 in the Ramond sector of left movers. The effect on Chan- 
Paton factors of open strings on D3 branes is given by the matrices 73 = i J2N C for the orbifold 
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action and CJ3 = l2N c for the orientifold action. The corresponding matrices for the D7 branes are 
77 = iJ2N ( and oj 7 = iJ 2 N r 

The massless spectrum of D3 branes involves a vector multiplet on the worldvolume A)i23 
and three complex scalar multiplets describing the transverse motion X45, Xqj, X^. Before the 
projection those describe the field content of M = 4 U (2N C ) super Yang-Mills, that in N = 2 
language involves a vector multiplet and a hypermultiplet in the adjoint representation. The 
orientifold action is 

A)123 — > -W3 ^0123 w 3 > 

X45 -> -UsX&Uz 1 , (17) 
-^67,89 — > <^3 -^67,89 w 3 1 • 

Therefore, the orientifold projection for the gauge field is 

V^A^KA^-Al), (18) 

so the projected gauge field is antisymmetric and spans an SO(2N c ) algebra. The field X45 is in an 
antisymmetric (adjoint) representation, while for the fields ^67,89 the orientifold action projects 
them to a symmetric representation. 
The Z2 action of the orbifold is 

^u;A)123 -> 73 £u;A)123 

V U X 45 -> 73^X 4 5 7 3 " 1 > (19) 

^^67,89 -> -73 ^0)^67,89 73 -1 - 

The transformations of ^0123 and X45 are identical and produce fields in the adjoint representation 
of U(N C ). The projection on ^67,89 produces fields in a two-index symmetric representation. More 
explicitly, for the gauge field the projection is 

V^VuAp = \ (V^ + J 2 N c P u A li J-^ e ) . (20) 

The resulting theory is a N = 2 U(N C ) theory with a symmetric hypermultiplet. If one considers 
the orientifold projection alone, the theory is projected to N = 2 SO(2N c ) super Yang-Mills with 
a hypermultiplet in the two-index representation, we can think of this theory as the analog of the 
SO(2N c ) gauge theory of the QCD case. 

The D3/D7 spectrum is initially described by two 2iV c x 27Vf chiral multiplets H A describing 
strings from D3 to D7 branes and the reversed strings Ha = eABH B ^ . The orientifold and orbifold 
actions are as follows 

H A -> -ieAB {u3H B u^)* , V„H A -> jsV^H^ 1 . (21) 

Therefore, the projections acting on flavor fields are 

V U H A = \{H A + e AB (H B )* J-^ f ) , VJ> U H A = \ (V^H A + J 2Nc V w H A J~^ . (22) 

The resulting massless field is a,N = 2 hypermultiplet in the (iV c , Nf) representation, or Nf flavors 
in the fundamental representation of the U (N c ) gauge group. In the theory obtained from the 
orientifold projection alone there are Nf hypermultiplets in the fundamental representation of the 
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S0(2N C ) gauge group. Although the maximal possible flavor group is U(2Nf), in the theory at 
hand it is actually reduced to USp(2N{), due to the coupling between the chiral components of 
the hypermultiplets with the chiral component of the vector multiplet in the superpotential 

W ~ HXH. (23) 

Since X is in the adjoint of SO{2N c ), flavor indices in the superpotential are contracted with an 
antisymmetric form, which is invariant under a USp{2N{) C U(2N{) subgroup. 

The massless spectrum of D7 branes, that describes the BPS sector of flavored operators, is 
split between vector fields in the 0123 and 6789 directions, A)i23 and -46789, and a scalar field in 
the 45 directions, X45. Transformations act as 

-4 123 -OJ7 A% 123 l , VloA 012 3 -> 77^^0123 77 _1 . 

X 45 -»• -u 7 XJ 5 lo^\ V U X& -»• 77 7^X 45 77" 1 , ( 24 ) 

-46789 -> -W7 ^6789 w 7 ) ^^6789 -> ~77 ^^6789 7f 1 • 

Since the 8d Poincare invariance is broken in the worldvolume of the D7 branes, the projection 



will be different for modes with dependence on the 6789 directions. The action (24) for ^4oi23 and 
X45 is valid for parity even modes while the action for ^6789 is valid for parity odd modes. This 
agrees with the A0123 and X45 components being scalar in the 6789 directions and ^789 being 
a vector component. The flavor group is J7(iVf), but the spectrum of BPS operators is different 
from the original theory since the hypermultiplet in the D3 sector is in the two-index symmetric 
representation and not in the adjoint. 

The holographic dual description is type IIB string theory on AdS§ x RP 5 , with D7 probe 
branes that sit on top of 07 planes wrapping a RP 3 C RP 5 cycle. The AdS§ x RP 5 geometry 
can be understood using a different basis of transformations. The 07 action is Q7 = QR^(— 1) Fl , 
while the Z2 singularity acts as a i?6789 reflection on the geometry. Since 07 planes and 03 planes 
have the same effect on Ramond forms (cf . [29] ) , the combined action is equivalent to the action of 
an 03 plane $^3 = i?6789^7 = ^i?456789( — l)^. The action of the 03 plane on AdS$ x S 5 is known 
to give the RP 5 geometry, since it acts as a reflection on the space transverse to the D3 branes 
[30J. From the T-dual perspective this geometry without the 07 orientifold can be constructed 
from a stack of D4 branes sitting on 04 _ or 04 + , giving holographic duals with orthogonal or 
symplectic gauge groups. 

3.2 From isospin to baryon chemical potential 



The dynamics of the probe D7 branes in the D3 background are determined by the DBI action, 

S DBI = -T 7 J d s £ Tr^-det (G + 2-ira'F), (25) 

where £ are the world- volume coordinates, G is the pull-back of the spacetime metric to the world 
volume and F is the field strength of the gauge fields on the brane. The 2iVf D7 branes in the 
U(2N C ) theory are wrapping an S ?J C <S 5 . Writing the AdS§ x S 5 metric as 

ds 2 = Via , dx^dx u + — £ , ( 26 ) 

i=4 
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the D7's are localized at yg = yg = and extend along all the other directions. As we have 
explained, the full projection identifies points in the geometry that map to each other under a 
reflection y l — > — y\ 

An isospin chemical potential in the U(2N C ) field theory is described by a background gauge 
field on the D7 brane 

A = i(j,J 2Nf . (27) 
More generally, this will be taken as the boundary condition for Aq . Notice that this configuration 



survives both projections (24). To highlight the effect of the projection let us write the gauge 



potential on the D7 brane as 

/i_ — / 

x Ct H' 

where H, H' and C are N{ x Nf Hermitian and general complex matrices, respectively. The 



Ao = ( 5 £ ) , (28) 



orientifold projection in (24) implies the conditions 

H' = -H T , C = C T , (29) 
so Aq is in the adjoint of a U Sp(2Nf) group, 

J 2Ni A + A$J 2N{ = 0. (30) 



The orbifold projection based on the transformation (24) imposes the conditions 

H = -H T , C = -C*. (31) 

So H is reduced to a purely imaginary antisymmetric matrix and C to a purely imaginary sym- 
metric matrix. The combination Aq = —iC — H belongs to the adjoint representation of a U (Nf) 
group. One can check this by doing a global transformation Aq — > UAqU^ with elements of the 
unbroken gauge group 

Uuj 7 U T = u> 7 , £777 U j = 77. (32) 



In terms of the U(N{) gauge field, the configuration (27) maps to 



Ao = nl Nt , (33) 

which corresponds to a baryon chemical potential. This shows that indeed the isospin chemical 
potential is projected to a baryon chemical potential. Notice that after the orientifold projection, 
because the gauge group SO(2N c ) is real, the isospin chemical potential is equivalent to the 
"baryon number" chemical potential. 

3.3 Validity of the equivalence 



The definition of the non-Abelian DBI action without derivatives is ambiguous, as it is possible 
to use [Dfi,D v ]F ai p = [F^ u ,F a/ 3] to convert between derivatives and field strengths. Nevertheless, 
a possible approach is to define the non-Abelian DBI action by giving an ordering prescription for 
the trace when the gradients of the field strength are small [38J. In this regime, the DBI action 
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Figure 3: Sketch of the phase diagram in the chemical potential [x and temperature T plane, for 
either baryon or isospin chemical potential. The phase with spontaneous isospin breaking only 
exists for the case of an isospin chemical potential. The large N c equivalence relating baryon and 
isospin chemical potential is valid in the confined and deconfined phases, but not in the region 
where isospin can be spontaneously broken. 



does not change under the projections. We show this more explicitly in Appendix [B] Under this 
assumption, we obtain the same results in the original and in the daughter theories, as long as 
all the components that are not invariant under the projection are zero. Notice that if restricted 
to questions about the ground state configuration with or without charge density, there are no 
ordering ambiguities; the configuration is Abelian so the separation between derivatives and field 
strengths is well-defined and any ordering prescription leads to the same results. 

From the perspective of the dual field theory, this means that the parent and the daughter 
theories are equivalent in the large- N c limit. In other words, the equivalence between the original 
and the daughter theory holds as long as the solution to the equations of motion obtained from 
the DBI action and the boundary conditions are invariant under the projection symmetry. We 
can fix the boundary conditions, but there is still the possibility that the correct solution breaks 
spontaneously one of the symmetries we have used to define the projection. This is a dynamical 
question, but fortunately one that can be answered in this context. 

The phase diagram of D7 branes in the D3 background geometry was studied in ref. [31] for 
a baryon chemical potential, and also in refs. \32\ [36]. [33] for an isospin chemical potential^] We 
have sketched the phase diagram in Figure [3] For non-zero quark mass m there is a phase with 
no charge density, where quarks can form bound states of mass ~ m/vA and there is a discrete 
spectrum of low-spin mesonic states. For large enough chemical potential and/or temperature 
there is a phase transition^] and a finite charge density (baryon or isospin) appears. In this phase 
the quarks do not form bound states and the spectrum is continuous. The curves of the phase 
transition in the (hb,i,T) plane are the same for both baryon and isospin chemical potential, as 

"Notice that we are considering a baryon chemical potential in a theory where two of the adjoint fields have 
been changed to two-index symmetric representations. However, both theories are equivalent in the large- N c limit. 
In the gravity side we have geometries with the same metric but different topologies, so solutions to the classical 
equations of motion are the same in both cases. 

12 The phase transition is first order at zero chemical potential, but changes to 3rd order at a tricritical point [34] 
at finite T and p. At zero temperature the transition becomes second order |35| . 
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we expect from the equivalence. The equivalence also holds in the finite density phasej^j but for 
an isospin chemical potential there is also a symmetry broken phase where vector mesons charged 
under the £7(1) isospin symmetry condense [36, 37] . The broken phase appears for values of the 
chemical potential larger than the scales given by the temperature and the mass of the mesons, 
which is proportional to the quark mass in these models. This implies in particular that at zero 
temperature and in the chiral limit the equivalence will fail. However, for massive quarks and 
small chemical potentials the equivalence will still be valid. 



3.4 Nf/N c corrections 



It is interesting to study the question of Nf/N c corrections in the theories with holographic 
duals and compare with QCD-like theories. Up to now, we have been working in a probe approxi- 
mation for the flavor branes. From the field theory perspective this is equivalent to working in the 
't Hooft limit of Nf/N c — > 0, with no quark loops in planar diagrams. To find the effect of Nf/N c 
corrections, one should compute the backreaction of the brane on the geometry. In principle, 
orientifold planes also produce a backreaction on the geometry of the same order; from the point 
of view of the field theory dual this corresponds to 1/N C corrections not associated with flavor. 
However, in the limit where the number of flavor branes is large, but still much smaller than the 
number of color branes, N c 3> Nf 3> 1, one can neglect the backreaction of the orientifold to 
leading order. From the field theory perspective, one can do a double expansion in Nf/N c and 
1/N C , where N{/N c corrections appear in planar diagrams with any number of flavor loops and 
1/N C corrections correspond to non-planar diagrams. 

The comparison of N{/N c corrections is a bit subtle. On the one hand, the DBI action is the 
same for both the normal configuration dual to isospin chemical potential, and the orientifolded 
configuration dual to baryon chemical potential. Wess-Zumino terms with an even number of field 
strengths are also the identical. However, terms with an odd number of field strengths can differ. 
Consider, for instance, the coupling of the D7 to the RR potential Cq, 

Swz = m [ C 6 Atr(F). (34) 
JD7 

Before introducing the orientifold, tr(i ? ) = 0, since the gauge field on the brane is non-Abelian 



(27). The projection does not affect to the value of the trace, but for a D7 brane with an Abelian 



configuration (33), tr(F) ^ in general, so this term is different in both cases and the backreaction 
of the D7 branes will also be different. Therefore, Nf/N c corrections in the parent and daughter 
theories will be different in general and the equivalence only holds in the strict 't Hooft limit. 



From the field theory perspective, we saw in section (2.3) that the condition (16) is not 
satisfied in the projection to a baryon chemical potential. For the supersymmetric theory we 
have to check two projections: both the orbifold projection, and the orientifold projection. The 



orbifold projection "P 7 mapping SO(2N c ) — > U(N C ), whose action is given eqs. (20) and (22), 
naively seems to be regular. However, the projection is written in a basis where fundamental 
fields are in a real representation, in particular it is valid for Majorana fermions, while in the 
daughter theory fermions will be Dirac. In terms of a complex representation with Dirac fermions 



In ref. (35] it was observed that there is a 4-fold symmetry in the (hb,Hi) plane. 
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one recovers the same kind of projection as in (|3j), which is not regular. Therefore, we also find 
from the field theory side that the equivalence will not hold for Nf/N c corrections. 

One can show indirectly that the orientifold projection Vu mapping U(2N C ) — > SO(2N c ) in 



eqs. ( |18| ) and (22) is regular. For this, notice that there is a regular projection from SO(AN c ) —> 
U (2N C ) with isospin chemical potential using a Z4 subgroup. This projection can be extended to 
SO(4N c ) — > SO(2N c ) using a D4 dihedral subgroup, with Z4 C -D4. Indeed, using the orbifold 
action -)• Jm c A^ to project 50(4iV c ) -)• U(2N C ) and 

^ - ( -L ) ■ <»> 

to further project to SO(2N c ) through the action — > L^^A^L^ , one generates a regular 



representation of D4 We can then group flavor fields in a 4N C x 4iVf real matrix 



_ , Hi + Hl H 2 + H% . 
1 i{Hl-H t ) i(H*-H 2 ) I ' 



and perform the orientifold projection ([22]) in this basis as 

T^ft = I (ft + £4iV c ft%4) » (37) 

where 

n - = ( 4 ) • (38) 

Since tr^jVf = and = l4j\r f , this shows that the orientifold projection in the flavor sector 
is also regular. Notice that the projection acting on fields on the flavor D7 branes is determined 
by tl^Nf , so from the perspective of the field theory on the D7 branes this is a regular projection. 
In the holographic duals to both the parent U(2N C ) and daughter SO(2N c ) theories the gauge 
field configuration on the D7 brane is traceless, so the issue of terms with an odd number of field 
strengths does not arise in this case. 



4 Conclusion 

We have used a holographic construction to demonstrate a large- iV c equivalence between 
theories with baryon and isospin chemical potential in the 't Hooft limit. The equivalence is valid 
in the region of the phase diagram where neither isospin nor baryon symmetry are spontaneously 
broken. The allowed region contains a small temperature/chemical potential phase with no charge 
density and mesonic bound states, and a large temperature/chemical potential phase with a finite 
charge density and a continuous spectrum. Although the charge density vanishes in the low 
temperature phase in the classical supergravity approximation, it will be nonzero when one takes 
Hawking radiation from the black hole into account. From the perspective of the dual field theory, 
this indicates that any charge density is suppressed in the large- iV c limit, compared to the charge 
density that is present in the high temperature phase. In other words, in the low temperature 
phase the equivalence we have presented relates the leading order large- iV c behavior of "vacuum" 

14 The elements of the group are {l4jv c , — l4iv c , J4jv c , — ^4jv c , i4iv c , — i4jv c , J4jv c l/4jv c , -£<4jv c J4jv c }. 
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thermodynamic properties of the theory with either baryon or isospin chemical potential, but 
it does not provide information about other details, such as properties of the thermal gas of 
mesons and baryons, which are 1/N C suppressed. The large iV c equivalence does relate non-trivial 
dependence on temperature and chemical potential in the deconfined phase (with unbroken isospin 
and baryon number), where there is temperature and chemical potential dependence at leading 
order in N c . 

Accepting the validity of gauge/string duality, our analysis suggests that large N c equivalences 
relating theories with differing chemical potentials may be valid more generally (when appropriate 
symmetry realizations hold). Exactly the same projections which we have used for supersymmetric 
theories can also be used to relate U(2N C ) QCD with 2Nf flavors and isospin chemical potential to 
U(N C ) QCD with Nf flavors and baryon chemical potential. However, because there is no known 
gravitational dual of QCD, different methods are needed to construct a purely field theoretic proof 
which applies to this case. We hope to revisit these points in the future. 
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A A quick introduction to the sign problem 

Consider pure Yang-Mills theory. In lattice Monte-Carlo simulations, one generates field 
configurations with probability weight e~ SyM /Zym (times Haar measure), where Sym is the 
Euclidean action and the partition function Zym = / dA^e~ SYM ^. Expectation values are 
approximated by taking the average over many configurations generated by a Markov chain: 

k 

(O) = lim \ V 0[A$] = — !— f dAu 0[A] e ~ SYM ^ . (39) 
fc^oo k f-f F Z Y m J 
i=\ 

Here k is the number of lattice field configurations generated in the simulation and i = 1, ■ ■ ■ , k is 
a label distinguishing them. In the sequence of configurations, more likely configurations appear 
more often; this is known as the "importance sampling" . 

The key assumption, which is valid for pure gauge theories, is that the the weight e~ SvM /Zym 
is real and positive, so that the integrand of the functional integral may be regarded as a probability 
measure. This condition can be broken when there are fermions in the theory. 

To deal with fermions in a lattice theory, one performs the Grassmann integral by hand. For 
example, the expectation value of the chiral condensate can be expressed as 

_ !dA,Trp-M-^p[A]-e-^ 
WW JdA^detp[A]-e-^M[A] ' ^ 

where If) denotes the lattice Dirac operator. If det p[A] is real and positive for any gauge field A^, 
then one may simulate this system by using the effective action S e ff[A] = 5yAf [i] — log det $>L4]. 
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However, the determinant det$)L4] can, in general, be complex, and then standard Monte-Carlo 
techniques cannot be applied. This is the so-called "sign problem" (or more properly "phase prob- 
lem"), the word "sign" referring to a possible negative sign of the determinant. QCD with baryon 
chemical potential suffers from the sign problem, while QCD with isospin chemical potential and 
the SO(2N c ) gauge theory with fermion chemical potential are "sign-free". 

One standard approach for dealing with the sign problem is the so-called "reweighing" method. 
Consider the phase quenched ensemble with weight | det Jp\e~ SyM . When the number of flavors 
is even, the phase-quenched version of QCD with baryon chemical potential is identical to QCD 
with an isospin chemical potential. If (• • • )b denotes expectations in QCD with a baryon chemical 
potential, and (•••)/ expectations with an isospin potential, then it is immediate that 

«°»-^ 

where e lv is the phase of the fermion determinant in the presence of a baryon chemical potential, 
e %ri = det If) /| det^)|. Because one can apply standard Monte-Carlo simulation techniques to the 
phase-quenched ensemble, one can in principle evaluate the expectation value in the full theory by 



computing both numerator and denominator in the identity (41); this is the reweighing method. 
However, in practice reweighing works only when the phase does not fluctuate violently. Phase 
fluctuations grow as the chemical potential \ib is increased, and as the lattice volume V grows. 
The logarithm of the fermion determinant is extensive, In det JJ) = 0(V), and (with a non-zero 
baryon chemical potential), so is its imaginary part, r/. This implies that both numerator and 
denominator of eq. ( |41[ ) vanish exponentially in the thermodynamic limit, making their estimation 
via Monte Carlo methods increasingly problematic as the volume V — > oo. 

B Projection of the DBI action 

The DBI action of a D{d— 1) brane may be expressed as 

S DBI = -T d ^ J d d £ S'Tr 2NiV / -det(G + 2ira'F), (42) 

where we use the notation S'Tr2jv f for a trace of 2iVf x 2iVf matrices with a predetermined ordering 
prescription that is, however, unknown. Up to F 4 terms, it should coincide with the symmetrized 
trace prescription [38] , 

For a d-dimensional brane, one can write the determinant as 

det (G + 2WF) = | e w-w e "i-^(G w „ 1 + 2-aa'F^) ■ ■ • (G w „ d + TnctF^) , (43) 

where G^ lUl is the pullback of the metric and is proportional to the identity matrix. The field 
strengths F^ lUl are proportional to the generators of the gauge group on the brane. Since this 
is a matrix product one has to define the order, we will not assume a particular ordering in the 
following. We can extract the metric factors as 

det (G + 2ira'F) 

= h e^-V^G, lQl ■ ■ ■ G„ dad (8% + 2™'G a ^F PlV1 ) ■ ■ ■ (8% + 2na'G a ^F^ d ) 
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= i (det G) 6" 1 '" w e ai ... Qd (^ + WG«^ m ) • • • (<5£ + 2WG^i^J . (44) 
Alternatively, one can write the determinant as 
det [G + 2ita'F) 

= ar e^-^e^-^G^^ • • • G adVd (6% + 2™'^^) ■ ■ • («5£ + Wff^^J 

= £ (det (?) e^-^e ai ... ad (^ + WG^i^J • • • (5 a / d + 2na/G a *f>* F^) . (45) 

Since G a "Fp v = —G a ^F u p, after relabeling indices v% /ij in (45) and comparing with (44), one 



finds that det [G + 27ra'F) = det (G — lira 1 F). Therefore, only even powers of the field strength 
F appear in the expansion of the determinant, which can be seen as a manifestation of charge 
conjugation invariance. Taking the square root and defining <jn (T N ) as a possible ordering of N 
generators T appearing in the tracePjwe have 



S"Tr 2 jv f V- det i G + 2vra'F) = V- det G 

( k [d/2]\ 

2N { + J2(a>r EE II £ ^ q n q -NC N t n2 ,.., nk tra 2N (F^F^...F 

N>1 <T2N k>\ \q=ln q =oJ 



2n k . 



(46) 



Where we have suppressed spacetime indices and denote F n as a product of n field strengths 
appearing in the determinant. The largest possible power is d if d is even or d — 1 if d is odd. 
The action has to be Hermitian, this implies that given some ordering 027V, the reversed ordering 
ctJn also appears with the same coefficients. For instance, if we have tr (F± F2 ■ ■ ■ F n —\ F n ), the 
Hermitian conjugate is tr (F n F ra _i • • • F2F1). 



Although expression (46) takes the form of an a' expansion, we are overlooking other a' correc- 
tions involving derivatives of the field strength and a' corrections that depend on the background 
metric. The former can formally be included by allowing DF factors inside the traces; many new 
terms (with unknown coefficients) would appear but the basic structure would not change if dis- 
crete symmetries (C,P,T) are not broken. We will just assume that gradients of the field strength 
are much smaller than their magnitude. Corrections to the geometry are more problematic, as 
a' corrections are, in general, different in the presence of an orientifold plane. But since we are 
both in the supergravity approximation and in a probe limit, we can neglect these corrections so 
the coefficients of the expansion are not affected. Under these assumptions the DBI action of the 
daughter theory is the naive projection of the original theory, as we will now show. 

Let us write the field strength of the gauge field on the D-brane as 



H C 
Ct H' 



(47) 



Here and in the following, we suppress spacetime indices. Under a global U (2N{) transformation 
taking F — > PFP^ 1 , with P defined in ([6]), the field strength transforms to 

PFP~i = l( H + H' + i(C^-C) H-H' + i(tf + C)\ 

2 V H-H' -i(C + Ct) H + H' — i(C$ — C) 1 



15 Since (jjv is defined for a trace, it maps to an element of the quotient of the permutation group over the cyclic 
group S n /Cn, which corresponds to the conjugacy classes of S„. These can be classified using Young tableaux. 
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After applying the conditions (29) and (31) which follow from the projection, the transformed 
field strength is block-diagonal, 



PFP- 



H-iC 






H + iC 



F T 




-F 



(49) 



Here F is the U(N{) gauge field of the daughter theory. Any power of the field strength F n has 
the same block-diagonal form after the global transformation, so for I 



L factors, the trace 



is 



tr °l n **" )= tr °l [ p n ptk p_i = tr m n^ T ) 4 + tr °m il-^ 



II^ 4 +(-l) i traJn^ 4 



(50) 



Since L is always even in the expansion (46), the phase factor is trivial (—1) = (—1) = 1. 
For each ordering a 2 N m (46) we get two terms, one corresponding to the same ordering for the 



U (iVf) gauge field and another one corresponding to the reversed order. From the Hermiticity of 
the action in the U(2N{) theory, we should have another contribution in the projected action that 
is exactly the same but whose origin is a term with reversed order cf^n- Adding the two together, 
we have a Hermitian action where the traces are projected to 



tr o 2N {F 2ni F 2na • • • F 2nk ) -> 2 tr a 2N (F 2ni F 2n2 • • • F 



Finally, using expression (46), 



S'Tr 2Nf V- det ( G + 2vrQ'F) -»• 2 S'Tr 



•jv f 



det ( G + 27ra'F 



(51) 



(52) 



This proves the equivalence of the ordering prescription in the original and the projected actions. 
Comparing with the action (42), the tension of the D-brane in the U(2N{) theory is half the 
tension of the U(N{) theory. However, this is compensated by the volume of the internal space, 
that is halved when we project from the S 5 to the RP 5 geometry. 
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